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1 .  Introduction 


This  is  the  final  report  on  the  research  program  on  "Radar  Waveform  Synthesis 
for  Target  Identification"  supported  by  the  Naval  Air  Systems  Command  under 
Contract  N00019-84-C-0190,  and  it  reports  the  progress  for  the  period  of 
July  31,  1984  to  August  31,  1985. 

The  purpose  of  this  research  is  to  develop  a  new  scheme  of  radar  discrimination 
and  identification.  The  new  scheme  is  based  on  the  natural  frequencies  of  the 
target.  It  consists  of  synthesizing  aspect-independent  discriminant  signals, 
called  Extinction-pulses  (E-pulses*)  and  single-mode  extraction  signals  which, 
when  convolved  numerically  with  the  late-time  transient  response  of  an  expected 
target,  lead  to  zero  or  single-mode  responses.  When  the  synthesized,  discriminant 
signals  for  an  expected  target  are  convolved  with  the  radar  return  from  a 
different  target,  the  resulting  signal  will  be  significantly  different  from  the 
expected  zero  or  single-mode  responses,  thus,  the  differing  targets  can  be 
di scriminated. 

The  complex  natural  resonant  frequencies  of  a  radar  target  are  aspect 
independent  features  of  its  transient  electromagnetic  response.  A  number  of 
researchers  have  recently  attempted  to  discriminate  among  various  targets  by 
extracting  those  natural  frequencies  from  late-time  transient  radar  returns. 

Since  extraction  of  natural  frequencies  from  late-time  target  responses  is 
an  inherently  ill-conditioned  numerical  procedure,  very  large  S-N  ratios  are 
required  in  the  transient  return.  It  has  therefore  been  concluded  that  this 
method  for  the  direct  discrimination  of  differing  target  is  impractical. 


* 

The  E-pulse  is  similar  to  the  K-pulse  studied  by  other  workers 


[1-2]. 
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Our  discrimination  scheme  differs  significantly.  Synthesis  of  the  discriminant 
signals  requires  only  knowledge  of  the  natural  frequencies  of  various  expected 
targets.  The  latter  natural  frequencies  are  measured  in  the  laboratory  where 
they  are  extracted  from  the  late-time  pulse  responses  of  target  scale  models. 

The  numerically  ill-conditioned  natural  frequency  extraction  procedure  need 
therefore  be  applied  only  to  target  responses  measured  in  a  controlled  S-N 
environment.  Synthesized  discriminant  signals  based  upon  those  laboratory 
measurements  are  stored  as  computer  data  files,  and  subsequently  convolved 
numerically  with  actual  transient  target  radar  returns.  Since  the  latter 
convolution  operation  is  numerically  well  conditioned  (a  smoothing  integral 
operator),  the  S-N  requirements  for  the  actual  radar  return  are  significantly 
relaxed. 

Another  observation  made  in  the  course  of  our  study  is  worth  noting. 

It  is  common  thinking  among  many  researchers  that  radar  detection  utilizing  the 
late-time  transient  radar  reutrn  may  not  be  practical  because  it  contains  little 
energy;  most  energy  is  associated  with  the  early-time  part  of  that  return.  This 
thinking  may  be  true  for  very  low-Q  targets.  Fortunately,  for  most  space 
vehicles,  such  as  rockets  and  aircrafts,  these  targets  are  not  exactly  low-Q 
structures.  There  is  sufficient  energy  contained  in  the  late-time  returns  of 
such  targets,  as  can  be  evidenced  from  our  measured  responses  of  complex  targets 
as  discussed  in  Section  3. 

Under  the  sponsorship  of  Naval  Air  Sy|tems  Command,  the  research  program 
has  progressed  steadily  over  the  past  few  years  and  so  far  we  have  a  good 
understanding  of  the  basic  principle  and  have  demonstrated  the  feasibility  and 
applicability  of  our  scheme  in  discriminating  between  complex  radar  targets 
[3-6].  It  appears  that  this  scheme  has  a  good  potential  to  be  a  useful  and 
practical  method  for  radar  detection  in  the  future. 


In  this  report  we  outline  the  progress  made  over  the  past  year.  In  Section  2, 


the  basic  theory  on  E-pulse  and  single-mode  extraction  signals  is  outlined  based 
on  the  time-domain  and  the  frequency-domain  analysis.  In  Section  3,  experimental 
results  of  the  scheme  when  applied  to  complex  targets  are  given.  In  Section  4, 
shaping  extraction  signals  by  proper  choice  of  basis  functions  is  discussed. 

In  Section  5,  new  methods  for  extracting  the  natural  frequencies  of  a  complex 
target  from  its  measured  pulse  response  are  presented.  The  future  plans  are 
given  in  Section  6.  Two  recent  papers  published  by  us  are  included  in  Appendices. 

2.  Theory  on  E-Pulses  and  Single-Mode  Extraction  Signals 

Basic  theory  on  E-pulses  and  single-mode  extraction  signals  is  briefly 
outlined  in  this  section. 

2.1.  Time  domain  analysis 

Assume  that  the  measured  time-domain  scattered  field  response  waveform  of  a 
conducting  radar  target  can  be  written  during  the  late-time  period  (t  >  T  )  as 
a  sum  of  damped  sinusoids 

N  at  /  t  \ 

r(t)  «  i  a  e  cosU  t  +  *  )  t  >  T^  K  1 

n»l 

where  aR  and  $n  are  the  aspect  dependent  amplitude  and  phase  of  the  n'th  mode, 
sn  *  °n  +  Jwn  the  asPect  Independent  natural  frequency  of  the  n’th  mode,  and 
only  N  modes  are  assumed  to  be  excited  by  the  incident  field  waveform.  Then, 
the  convolution  of  an  E-pulse  waveform  e(t)  with  the  measured  response  waveform 
becomes 

rTe 

c(t)  *  e(t)*r(t)  *  I  e(t')r(t  -  t ’ ) dt ' 

J  0 

N  o  t 

*  t  a  e  n  [Acos(u  t  +  *  )  +  B  s1n(u>  t  +  ♦  )]  (2) 

,  n  n  n  n  n  n  n 


for  t  >  T,  *  T  ♦  T 


where 


f e  -at'  ^ 

An  *  I  e(t'  )e  cosu^t'dt' 

V  (3) 

rTe  -0  t' 

B  *  e(t')e  sinuj  t'dt' 

n  Jo  n  J 

and  Tg  is  the  finite  duration  of  e(t). 

Two  interesting  waveforms  are  now  considered.  Constructing  e(t)  to  result 
in  c(t)  *  0,  t  >  Tl,  requires 


An  ~  B  =  0 
n  n 


1  <  n  <  N 


In  addition,  e(t)  can  also  be  constructed  so  that  c(t)  is  composed  of  just  a 
single  mode.  In  this  case  e(t)  is  termed  a  "single  mode  extraction  waveform". 
If  the  phase  of  c(t)  is  unimportant,  e(t)  can  be  constructed  by  demanding 


A  =  B  *  0 
n  n 


1  <  n  <  N,  n  ^  m 


to  excite  the  m’th  natural  mode.  On  the  other  hand,  requiring 


An  "  Br 


A  -  0 
m 


results  in 


1  <  n  <  N,  n  f  m 


c(t)  «  a  e  B  s1n(ui  t  ♦  *_) 
m  m  m  m 


and  requiring 


A  *  B 
n  n 


1  <  n  <-  N,  n  f  m 


B  =  0 


fc'TS 


The  E-pulse  resulting  from  (5)  is  termed  a  ”sin/cos“  single  mode  extraction 

waveform,  since  it  excites  both  sine  and  cosine  components  in  c(t),  while  (6) 

results  in  a  "sine"  and  (8)  in  a  "cosine"  single  mode  extraction  waveform.  With 

the  proper  normal izations  of  e(t)  (giving  A  *  B  ),  the  convolved  waveforms 

(7)  and  (9)  can  be  combined  to  yield  the  frequency  of  the  m'th  mode,  s  *  a  +  ju 

m  m 

2.2.  Frequency  domain  analysis 

The  convolution  of  the  E-pulse  waveform  with  the  measured  response  waveform 
can  also  be  written  in  the  form 
N  at 

c ( t )  =  i  an|E(sn)',e  n  cos^t  +  *n)  t  >  TL  (10) 


where 

E(s)  *  £*{e(t) } 


e(t)e’stdt 


is  the  Laplace  transform  of  the  E-pulse  waveform,  and 


(ID 


tii  *  <t  ♦  tan 
n  n 

where 

E1n  =  In,lE‘sn»  Ern  *  RetE<sn)l 
Now,  c(t)  «  0  for  t  >  T^  requires 


Ein 


(12) 


(13) 


In  rn  u 
or,  equivalently, 
E(sn)  -  E(s*)  -  0 


In  addition,  a  sln/cos  si 
E(sn)  •  E(s*)  *  0 


1  <  n  <  N  (14) 


1  <  n  <  N  (15) 

ngle  mode  extraction  waveform  can  be  constructed  via 
1  <  n  <  N,  n  i1  m  (16) 


while  a  sine  single  mode  extraction  waveform  requires 


It  is  easily  shown  that  the  frequency  domain  and  time  domain  requirements 
for  synthesizing  an  E-pulse  are  identical.  By  expanding  the  exponential  in 
(11),  one  can  show  that  (16),  (17),  and  (18)  are  equivalent  to  (5),  (6),  and 
(8) ,  respectively. 

One  benefit  of  using  a  frequency  domain  approach  comes  via  the  increased 
intuition  allowed  by  equation  (10).  When  an  E-pulse  waveform  is  convolved  with 
the  measured  response  of  an  unexpected  target,  the  amplitudes  of  the  resulting 
natural  mode  components  are  determined  by  evaluating  the  magnitude  of  the 
spectrum  of  e(t)  at  the  natural  frequencies  of  the  target  (a  result  of  the 
Cauchy  residue  theorem).  Thus,  the  E-pulse  spectrum  becomes  the  key  tool  in 
predicting  the  success  of  E-pulse  discrimination. 

2.3.  E-pulse  synthesis  (using  Frequency  Domain  Analysis) 

To  implement  the  E-pulse  requirements  it  becomes  necessary  to  represent  the 
waveform  mathematically.  Let  e ( t )  be  composed  of  two  components 

e(t)  =  ef(t)  +  ee(t)  (19) 


Here  e  (t)  is  a  forcing  component  which  excites  the  target,  and  e  (t)  is  an 
extinction  component  which  extinguishes  the  response  due  to  ef(t).  The  forcing 
component  Is  a  free  choice,  while  the  extinction  component  is  determined  by 
first  expanding  in  a  set  of  basis  functions 
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(20) 


where 


Fm(s)  3*{fm(t)} 
Ef(s)  (t) } 


(22) 


and  M  *  2N  Is  chosen  to  make  the  matrix  square. 

Two  types  of  E-pulses  are  now  easily  identified.  When  e^(t)  f  0,  the  forcing 
vector  on  the  right  hand  side  of  (21)  is  nonzero,  and  solutions  for  the  basis 
function  amplitudes  exist  for  any  choice  of  E-pulse  duration,  T  ,  which  does 
not  cause  the  matrix  to  be  singular.  In  contrast,  when  eF(t)  =  0  the  matrix 
equation  becomes  homogeneous,  and  solutions  for  e  (t)  exist  only  for  specific 
durations  Tg  which  are  calculated  by  solving  for  zeros  of  the  determinantal 
equation.  The  former  type  of  E-pulse  is  termed  "forced"  and  the  latter 
"natural".  Since  a  natural  E-pulse  has  no  forcing  component,  it  is  viewed  as 
extinguishing  its  own  excited  field. 
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A  very  useful  application  of  the  frequency  domain  approach  results  from  using 


pulses  as  the  basis  functions  in  (20).  Let 


fm(t> 


g(t  -  [m  -  l]u) 


(m  -  1  )i  <  t  <  mi 


el sewhere 


where  g(t)  is  an  arbitrary  function,  and  a  is  the  pulse  width.  Then 


F  (S) 


g(t  -  [m  -  l]^)e'SIdt 


*  F1 (s)eSAe'smA 


and  the  matrix  equation  (21)  can  be  written  for  the  case  of  the  natural  E-pulse 


1  Z,  ZZ  ...  zf-1  1  [a, 


1  ^  z; 


i  z;  (zy)‘ 


1  zj  <z*>‘  .. 


•  <Z5> 


,,2N-1 


where 


Zn  se 


Equation  (25)  is  homogeneous,  and  thus  has  solutions  only  when  the 


determinant  of  the  matrix  is  zero.  As  the  determinant  is  of  the  Vandermonde 


type,  the  condition  for  a  singular  matrix  can  be  calculated  easily  as 


P  *  1,2,3,. 


1  <  k  <  N 


<*  Z-  /  .*  .•  .-  .•  .•  J.  _.  . 


’hus,  the  duration  of  natural  E-pulse  depends  merely  on  tie  imaginary  part  of  one 
of  tne  natural  frequencies.  with  p  determined,  the  basis  function  amplitudes 
can  be  calculated  using  Cramer's  rule  and  the  theory  of  determinants  as 


-  <  1  Q 

am  ‘  '  ;  (2N  -  1)  -  (m  -  1', 


wn ere  R.  ^  is  the  sum  of  the  products  n  -  i  at  a  time,  without  repetitions,  of 
the  quantities  Z-j  ,  ZT,  Z?,...,Z£. 

Note  that  g(t)  does  not  appear  in  this  analysis,  and  thus  the  resulting  pulse 
amplitudes  are  independent  of  tne  individual  pulse  shapes.  However,  when 
di scriminating  between  different  targets,  g(t)  manifests  quite  importantly 
through  the  term  (s) . 

The  synthesis  of  single-mode  extraction  signals  can  be  carried  out  similarly 
based  on  (16),  (17)  or  (18).  The  only  difference  from  the  E-pulse  synthesis 
is  that  the  number  of  basis  functions  M  should  be  chosen  to  match  the  number 
of  equations  presented  by  (16),  (17),  or  (18). 

It  is  also  noted  that  the  E-pulse  synthesis  based  on  the  time  domain  analysis 
has  been  published  recently  [3]  (see  Appendix  1). 

3.  Experimental  Results  on  Complex  Radar  Targets 

In  the  preceding  section,  the  E-pulse  and  single-mode  extraction  signals  were 
synthesized  based  on  the  prior  knowledge  of  target's  natural  frequencies.  In  the 
case  of  complex  radar  targets,  this  information  is  difficult  to  obtain,  and  the 
synthesis  of  the  discriminant  signals  can  only  be  carried  out  by  a  combined 
experimental  and  theoretical  technique. 

Over  the  past  few  years,  we  have  developed  the  techniques  (which  can  be 
refined  further)  for  synthesizing  the  E-pulses  and  single-mode  extraction 
signals  for  complex  targets  based  on  measured  pulse  responses  of  their  scale 
models.  Through  the  convolution  of  these  synthesized  discriminant  signals  with  the 
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measured  radar  responses  of  the  targets,  we  have  definitely  demonstrated  the 
capability  of  our  scheme  to  discriminate  between  complex  targets. 

The  following  steps  are  used  in  synthesizing  the  discriminant  signals  for  the 
complex  targets: 

(1)  Measure  the  pulse  response  of  the  scale  model  of  the  target  at  various 
aspect  angles. 

(2)  Use  the  Fast  Fourier  Transform  to  obtain  approximate  values  of  natural 
frequencies  from  the  measured  pulse  response. 

(3)  Employ  the  continuation  method  [7  ]  and  the  natural  frequencies  obtained  from 
FFT  as  the  initial  guesses  to  calculate  accurate  values  of  natural  frequencies 
of  the  target.  At  each  step  of  the  algorithm,  the  condition  number  of  the 
regularized  problem  is  checked. 

(4)  With  the  natural  frequencies  of  the  target  determined,  the  theoretical 
technique  described  in  Section  2  is  then  applied  to  synthesize  the 
discriminant  signals. 

After  the  discriminant  signals  of  a  complex  target  were  synthesized,  they  were 
convolved  with  the  measured  radar  responses  of  the  target  at  various  aspect 
angles  to  check  the  workability  of  the  scheme. 

For  the  purpose  of  demonstration, two  complex  targets,  scale  models  of 
McDonnel -Douglas  F-18  airplane  and  Boeing  707  airplane  which  have  similar  sizes 
but  different  geometries,  are  used  here. 

Figures  1  and  2  show  the  measured  pulse  response  of  the  707  and  the  F-17 
scale  models,  respectively.  Each  model  is  constructed  of  aluminum  and  has  a 
geometry  as  indicated  in  the  figures.  Also  shown  are  the  dominant  natural 
frequencies  extracted  from  the  late  time  portion  of  the  response  using  the 
continuation  method.  These  frequencies  can  then  be  used  to  construct  natural 
E-pulse  and  single-mode  extraction  signals. 
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Pulse  basis  function  natural  E-pulses  of  minimum  duration  considered  for  each 
of  the  two  targets  are  shown  in  Figs.  3  and  d.  The  cosine  first-mode  and  the 
sine  first-mode  extraction  signals  for  the  707  airplane  are  shown  in  Fig.  5. 

Other  single-mode  extraction  signals  for  these  two  airplanes  are  not  shown 
here  for  brevity. 

Synthesized  E-pulses  and  single-mode  extraction  signals  for  these  two 
airplanes  are  then  convolved  with  their  measured  pulse  responses  shown  in 
Figs.  1  and  2.  The  convolved  results  are  depicted  in  the  following  figures. 

Figure  6  shows  the  convolution  of  the  E-pulse  for  the  F-18  plane  with  the 
measured  radar  response  of  the  F-18  plane.  The  convolved  output  shows  a 
strong  early-time  response  followed  by  a  "extingui shed"  (almost  zero)  late- 
time  response  as  expected.  Figure  7  shows  the  convolved  output  of  the  E-pulse 
for  the  707  plane  with  the  measured  radar  response  of  the  F-18  plane.  This 
convolved  output  shows  a  significant,  "unextinguished"  late-time  response 
implying  that  the  707  E-pulse  was  convolved  with  the  radar  response  of  a  wrong 
target  other  than  the  707  plane.  Figure  8  shows  the  convolution  of  the  E- 
pulse  for  the  707  plane  with  the  measured  radar  response  of  the  707  plane. 

The  convolved  output  shows  a  nearly  "extinguished"  late-time  response  because 
It  Is  the  right  target  for  that  E-pulse.  Figure  9  Indicates  the  convolved 
result  of  the  E-pulse  for  the  F-18  plane  with  the  measured  radar  response  of 
the  707  plane.  As  expected,  the  convolved  result  shows  a  significant 
"unextinguished"  late-time  response.  The  results  of  Figs.  6  to  9  clearly  show 
the  capability  of  target  discrimination  using  the  E-pulses  of  the  targets. 

To  enhance  the  certainty  of  the  target  discrimination,  the  single-mode 
extraction  signals  of  the  target  were  used  to  convolve  with  the  measured  radar 
responses  of  the  right  target  and  the  wrong  targets.  Figure  10  shows  the 
convolved  results  of  the  measured  radar  response  of  the  707  plane  with  the 
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first-mode  extraction  signal  for  the  707  plane.  The  left  figure  of  Fig.  10  shows 
the  angular  frequency  line  extracted  from  the  complex  convolved  output.  S’ nee 
this  line  is  almost  parallel  to  the  u,t  line  (of  the  707)  in  the  late-time 
period,  it  implies  that  the  measured  radar  response  is  from  the  right  target  of 
the  707  plane.  The  right  figure  of  Fig.  10  shows  the  damping  coefficient  line 
extracted  from  the  complex  convolved  output.  This  line  is  closely  parallel  with 
the  ^t  line  (of  the  707)  in  the  late-time  period,  implying  the  target  is  the 
707  plane.  Figure  11  shows  the  convolved  results  of  the  measured  radar  response 
of  the  707  plane  with  the  fourth-mode  extraction  signal  for  the  707  plane.  In 
this  figure,  the  extracted  angular  frequency  line  is  almost  parallel  to  the 
t  line  i  of  the  707  )  and  the  extracted  damping  coefficient  line  closely 
parallel  to  the  :^t  line  (of  the  707)  in  the  late-time  period.  This  implies 
that  the  radar  response  belongs  to  the  right  target  of  the  707  plane.  When  the 
measured  response  of  the  F-18  plane  is  convolved  with  the  first-mode  extraction 
signal  of  the  707  plane,  the  results  are  shown  in  Fig.  12.  The  extracted  angular 
frequency  line  is  not  parallel  to  the  ^t  line  (of  the  707)  and  the  extracted 
damping  coefficient  line  deviates  from  the  u^t  line  (of  the  707)  in  the  late¬ 
time  period.  This  indicates  that  the  measured  radar  response  comes  from  a 
wrong  target  other  than  the  707  plane.  Figure  13  shows  the  convolved  results  of 
the  measured  radar  response  of  the  F-18  plane  with  the  fourth-mode  extraction 
signal  of  the  707  plane.  The  extracted  angular  frequency  line  and  the  damping 
coefficient  line  deviate  from  the  .^t  line  (of  the  707)  and  the  n^t  line  (of  the 
707/,  respectively ,  in  the  late-time  period,  implying  that  the  radar  response 
belongs  to  a  wrong  target  other  than  the  707  plane. 

Results  depicted  in  Figs.  6  tc  13  definitely  confirm  the  capability  of  the 
target  discrimination  provided  by  the  method  based  on  the  concept  of  the  E-pulses 
and  single-mode  extractor  signals. 
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It  is  important  to  demonstrate  that  the  E-pulse  and  single-mode  extraction 
signals  of  a  target  are  aspect-independent .  The  following  examples  are  given 
for  this  purpose. 

Fig.  14  shows  the  convolved  outputs  of  the  E-pulse  of  F - 1 8  model  with  the 
pulse  responses  of  F-18  model  measured  at  (a)  0°  aspect  angle,  (b)  30°  aspect 
angle,  (c)  45°  aspect  angle  and  (d)  60°  aspect  angle.  In  each  case  the 
convolved  output  gives  an  insignificant  response  in  the  late-time  periods.  This 
implies  that  natural  modes  of  the  target  were  extinguished  by  its  E-pulses  and 
the  pulse  responses  used  in  the  convolution  belonged  to  the  F-18  model. 

Figure  15  shows  the  convolved  outputs  of  the  E-pulse  of  F-18  model  with  the 
pulse  responses  of  B 707  model  measured  at  (a)  0°  aspect  angle,  (b)  30°  aspect 
angle,  (c)  90°  aspect  angle  and  (d)  180°  aspect  angle.  It  is  observed  that  the 
convolved  outputs  give  large  late-time  responses  in  all  four  aspect  angles. 

This  indicate  that  the  F-18  E-pulse  was  convolved  with  the  pulse  response  of 
a  different  target. 

Figure  16  shows  the  convolved  outputs  of  the  E-pulse  of  B707  model  with  the 
pulse  responses  of  B707  model  measured  at  (a)  0°  aspect  angle,  (b)  30°  aspect 
angle,  (c)  90°  aspect  angle  and  (d)  180°  aspect  angle.  It  is  observed  that  the 
convolved  outputs  give  small  late-time  responses  in  all  four  aspect  angles  as 
expected.  This  implies  that  B707  E-pulse  was  convolved  with  the  pulse  responses 
of  the  same  target. 

Figure  17  shows  the  convolved  outputs  of  B707  E-pulse  with  the  pulse  responses 
of  F-18  model  measured  at  (a)  0°  aspect  angle,  (b)  45°  aspect  angle,  (c)  90° 
aspect  angle  and  (d)  180°  aspect  angle.  Large  late-time  responses  are  obtained 
in  the  convolved  outputs  for  all  four  aspect  angles,  signifying  that  the  pulse 
responses  used  in  the  convolution  belonged  to  a  target  different  from  B707  model. 
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Fig. 14.  Convolved  outputs  of  F-18  E-pulse  with  pulse  responses  of  F-18  model 
measured  at  (a)  0°  aspect  angle,  (b)  30°  aspect  angle,  (c)  45 
aspect  angle  and  (d)  60°  aspect  angle.  In  each  case,  a  small  response 
was  obtained  in  the  late-time  period. 
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Fig.  17.  Convolved  outputs  of  B707  E-pulse  with  pulse  responses  of  F-18  model 
measured  at  (a)  0°  aspect  angle,  (b)  45°  aspect  angle,  (c)  90°  aspect 
angle  and  (d)  180  aspect  angle.  In  each  case,  a  large  response  was 
obtained  in  the  late-time  period. 


The  results  of  Figs.  14-17  confirm  the  aspect  independency  of  the  E-pulse 
of  a  complex  target  and  the  feasibility  of  its  use  in  the  target  discrimination. 

4.  Shaping  Extraction  Signals  by  Proper  Choice  of  Basis  Functions 

In  the  synthesis  of  a  single-mode  extraction  signal,  it  is  important  for 
the  spectrum  of  a  single-mode  extraction  waveform  to  have  a  large  magnitude  at 
the  frequency  of  the  mode  to  be  extracted,  in  comparison  to  that  at  the 
frequencies  to  be  eliminated.  This  is  due  to  the  practical  problem  of  the  presence 
of  errors  in  the  measured  natural  frequencies  used  to  construct  the  single  mode 
extraction  signals.  Because  of  this  error,  the  modes  which  should  be  eliminated 
will  actually  be  extracted.  The  amplitudes  of  these  modes  are,  of  course, 
determined  by  the  amount  of  energy  in  the  extraction  signal  spectrum  at  the  modal 
frequencies.  Now,  if  there  is  a  relatively  small  amount  of  energy  at  the 
frequency  to  be  extracted,  these  extraneous  modes  may  make  a  relatively  large 
contribution  to  the  late-time  portion  of  the  convolved  response,  so  that  the 
expected  target  might  not  be  properly  identified.  This  is  especially  trouble¬ 
some  for  the  case  of  higher  mode  extraction,  where  the  smaller  damping 
coefficients  of  the  lower  order  modes  to  be  eliminated  allow  them  to  dominate 
the  expected  higher  order  mode  in  the  latter  part  of  the  late-time  region. 

Assume  that  we  wish  to  extract  a  single  mode  with  frequency  s  =  j  +  j*. 

Consider  using  the  damped  sinusoid  functions  for  the  basis  functions, 

Vl>  •  e  n  cos('mt  ♦  ;„)  (29) 

with  the  transforms 

t  sinh(1:  -  s>  -f 

(30) 


where  sm  *  3_  ♦  jZ  .  The  choice  of  s_  and  7  determines  the  usefu’ness  of  this 

set  of  basis  functions.  Alone,  each  F  (s'  has  a  peak  a’ong  a  i'ne  Re  s  *  constant 

near  *,  *  7  .  Thus,  choosing  Z  *  t  for  each  m  should  result  in  an  E-ouise  w’tn  a 
in  m 

spectrum  peaked  near  C. 

As  a  simple  numerical  example.  Fig.  18  shows  the  reconstruct 'on  z*  tne 
measured  pulse  response  of  a  McDonnel  Douglas  r-lS  aircraft  mode'  lee  : • ;  1 

The  three  dominant  modes  are  used  in  the  reconstructs  s.  1  [- 1  264  *  •  1  1?]«  1' 

$2  =[-0.126  j  7 . 32]  x  10^,  and  s^  =[-0.44'  *  ;  9 . 39)  *  '1*  Added  tc  t—s 
reconstructed  response  is  10  dB  of  random  noise  'he  noise  :ar  oe  .  e*ed  as 

perturbing  the  natural  frequencies  in  tne  '•espouse.  r-gure  '9  snows  t«c  *'r  et 

sine/cos  waveforms  constructed  to  excite  the  second  mode  z*  tne 
response.  The  first  waveform  is  synthesized  using  stangu'a'-  se  tas-.  ' 

with  a  forcing  function  chosen  to  be  another  -•ectangu ’ dr  :^'se  *•>  se-  r:  -a.e 

form  is  constructed  from  damped  sinusoi  d  funct'ons  «’t*  t"e  ‘  1  - : 

phases  given  in  the  figure.  Here,  the  forcing  ^unction  -s  '’,cser  1 
damped  sinusoid. 

Figure  20  shows  the  spectra  (magnitude!  0^  the  two  e«tra.. f  --  i-<ra 
plotted  at  Re(s)  x  While  the  spectrum  of  the  pulse  function  f  -:u  sr 

seen  to  have  its  signal  energy  concentrated  near  „  =  0,  the  spectrum  z* 
damped  sinusoid  based  E-pulse  has  a  majority  of  its  signal  energy  near  :ne 
frequency  to  be  excited.  The  result  of  this  energy  concentration  can  pe  seer 
in  Fig.  21,  which  shows  the  convolutions  of  the  two  E-pulses  with  the  noiseless 
reconstructed  response.  As  expected,  each  results  in  a  single  damped  sinusoid 
in  the  late-time,  but  the  late-time  response  due  to  the  damped  sinusoid  based 
E-pulse  is  seen  to  be  relatively  much  larger  Wore  importantly.  Fig.  22 
shows  the  result  of  convolving  each  E-pulse  with  the  noisy  response  of  Fig.  18. 

It  is  easily  seen  that  the  late-t’me  convoked  resrcrse  *rnr  ‘np  la-ipec 
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Reconstruction  of  measured  late-time  scattered  field  response  of 
McDonnel  Douglas  F-18  aircraft  model,  using  first  three  natural 
modes,  with  10  d8  of  random  noise  added. 
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Fig.  19.  Pulse  function  based  and  damped  sinusoid  based  forced  sin/cos 
single-mode  E-pulses  synthesized  to  eliminate  first  and  third 
modes  of  reconstructed  response,  while  exciting  the  second. 


relative  magnitude  of  spectrum  |E(s)| 


Convolution  of  pulse  function  and  damped  sinusoid  based  single 
mode  E-pulses  with  noisy  reconstructed  response. 


sinusoid  based  E-pulse  has  been  affected  much  less  by  the  perturbing  of  the 
natural  frequencies  of  the  reconstructed  target  response. 

There  is  a  need  to  study  the  uniqueness  of  E-pulses  and  the  possiblity  of 
minimizing  the  duration  of  the  E-pulse.  The  frequency  spectrum  of  an  E-pulse 
may  be  controlled  by  a  proper  choice  of  basis  functions.  Further  study  on  the 
E-pulse  may  make  it  more  insensitive  to  noise.  These  problems  will  receive 
attention  in  the  future. 

5.  New  Methods  for  Extracting  the  Natural  Frequencies  of  a  Complex  Target 
from  its  Measured  Pulse  Response 

In  our  scheme,  we  need  to  find  the  natural  frequencies  of  a  complex  target 
from  its  measured  pulse  response  before  the  E-pulse  and  single-mode  extraction 
signals  of  the  target  can  be  synthesized.  Since  the  classical  Prony’s  method 
is  known  to  be  extremely  sensitive  to  noise,  it  is  important  to  develop  other 
methods  for  this  purpose.  The  first  technique  developed  by  us  is  the  Continuation 
Method  which  has  been  published  [7]  (see  Appendix  2)  and  is  being  used.  This 
method  is  not  optimal  yet  and  we  plan  to  develop  other  different  methods.  Some 
preliminary  results  on  our  three  new  methods  are  given  below. 

The  measured  scattered  field  or  surface  current  response,  r(t)  of  a  target 
to  a  transient  pulse  waveform  is  assumed  to  be  composed  of  a  finite  number  of 
natural  oscillation  modes  in  the  late-time  period,  t  >  T  .  The  following 
three  techniques  attempt  to  extract  the  natural  frequencies  from  the  sampled 
representation  of  r(t). 

A.  E-Pulse  Method 

The  convolution  of  an  excitation  signal  with  the  measured  response  can  be 
written  as 

q(t)  *  e(t)  *  r(t)  (31) 
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The  natural  frequencies  contained  in  r(t)  can  be  extracted  by  constructing  e(t) 
as  an  E-pulse.  Then,  if  r(t)  is  indeed  a  sun  of  natural  modes,  q(t)  will  be 
identically  zero  in  the  late-time  of  the  convolved  response,  t  >  T^.  To  the 
extent  that  r(t)  will  be  contaminated  by  various  types  of  noise,  an  estimate 
of  the  natural  frequencies  can  be  obtained  by  minimizing  the  norm  of  the  late¬ 
time  convolved  response 

.  q(t)  2  =  :  [q(t.)]2  (32) 

i-1  1 

with  respect  to  the  parameters  0^...,  ^  used  to  construct  the  E-pulse  waveform, 
where  M  discrete  points  are  chosen  at  which  to  evaluate  the  norm.  The  natural 
frequencies  contained  in  r(t)  are  taken  to  be  ^  at  the  minimum  point. 

The  benefit  of  the  E-pulse  method  over  the  continuation  method  is  that  no 
initial  guesses  are  needed  for  the  phases  and  amplitudes  of  the  modes,  since 
they  are  not  involved  in  constructing  an  E-pulse.  This  also  means  that  the 
number  of  parameters  involved  in  the  minimization  is  reduced  by  1/2. 

B.  Discrete  Late-time  Minimization 

Rather  than  performing  a  minimization  with  respect  to  the  complex  frequencies 
used  to  construct  e(t),  this  method  performs  the  minimization  of  (32)  with  respect 
to  the  amplitudes  of  the  basis  functions  comprising  e(t).  Extraction  of  the 
natural  frequencies  in  r(t)  then  becomes  a  two  step  process.  The  minimization 
first  provides  an  estimate  of  the  E-pulse  waveform,  and  then  it  becomes  necessary 
to  determine  the  complex  frequencies  that  the  E-pulse  eliminates.  This  is 
accomplished  by  locating  the  zeroes  of  the  E-pulse  spectrum. 

If  a  rectangular  pulse  function  basis  set  is  used  to  construct  the  E-pulse, 
the  location  of  the  zeroes  of  the  spectrum  is  quite  simple.  The  spectrum  of  the 
E-pulse  becomes 
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where  F^(s)  is  the  Laplace  transform  of  the  first  pulse  function,  2  is  the 
pulse  width,  and  an  is  the  amplitude  of  the  n'th  pulse  basis  function.  Setting 
(33)  equal  to  zero  gives 
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which  is  just  a  polynomial  equation,  the  roots  of  which  give  the  complex 
frequencies  eliminated  by  the  E-pulse.  These  are  then  taken  as  an  estimate 
of  the  natural  frequencies  contained  in  r(t). 

The  benefit  of  this  method  over  the  continuation  method  is  the  same  as 
that  of  the  E-pulse  method  --  half  the  number  of  parameters  are  involved  in  the 
minimization,  and  there  is  no  need  for  initial  guesses  of  modal  amplitude  and 
phase.  An  initial  guess  for  the  basis  function  amplitudes  is  provided  by 
choosing  initial  guesses  for  the  natural  frequencies  and  constructing  an 
E-pulse  from  them.  However,  this  is  the  only  time  at  which  an  E-pulse  needs 
to  be  constructed,  which  should  represent  a  considerable  savings  in  execution 
time  over  the  E-pulse  method. 

C.  Moment  Method  Approach 

This  technique  is  also  a  two  step  process,  first  sovling  for  an  E-pulse 
using  the  method  of  moments  to  solve  the  integral  equation  created  by  setting 
the  late-time  convolved  response  to  zero,  and  then  extracting  the  frequencies 
from  the  E-pulse. 

Writing  the  convolved  response  (31)  in  integral  form,  and  assuming  that 
e(t)  is  an  E-pulse  results  in 

fT* 

q(t)  =  e(t')r(t  -  t')dt'  =0  t  >  Tr  (35) 

J  0  E 
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where  Tg  is  the  E-pulse  duration.  The  waveform  which  solves  this  integral 
equation  can  be  calculated  by  using  the  moment  method.  First  expand  e(t)  in  a 
set  of  basis  functions 
L 

e(t)  «  i  a  f  (t)  (36) 

z=l  1  1 

'hen  rather  than  forcing  (35)  to  be  satisfied  at  all  time  t  >  T^,  it  is 
multiplied  by  a  set  of  M  weighting  functions  Wm ( t )  and  the  moments  are  taken 


wnere  the  angle  bracket  represents  the  usual  inner  product. 

The  integral  equation  has  now  been  reduced  to  a  homogeneous  matrix  equation, 
where  M  is  chosen  to  be  equal  to  the  number  of  basis  functions  used  to  construct 
e(t).  A  solution  to  this  equation  demands  that  the  determinant  of  the 
coefficient  matrix  be  zero.  This  results  in  discrete  solutions  for  the  E-pulse 
duration  T  .  Although  a  search  is  required  for  these  durations,  the  amount  of 
complexity  Involved  is  much  less  than  that  of  either  the  E-pulse  method  or  the 
direct  late-time  minimization  method,  since  now  only  one  parameter  is  involved. 

Once  the  E-pulse  waveform  is  determined,  the  complex  frequencies  eliminated 
by  the  E-pulse  are  taken  as  estimates  of  the  natural  frequencies  in  r(t).  If 
e(t)  is  constructed  using  pulse  basis  functions,  then  the  detrmination  of 
these  complex  frequencies  proceeds  exactly  as  in  the  direct  late-time 
minimization  method,  by  the  solution  to  a  polynomial  equation. 

Two  choices  of  weighting  function  Wm(t)  used  in  (37)  have  been  investiated. 
Using  delta  functions  is  equivalent  to  point  matching,  and  requires  that 
equation  (35)  hold  at  discrete  points  in  the  late-time.  Using  rectangular 
oulse  functions  requires  that  an  average  value  of  the  late-time  convolved 
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response  be  zero,  and  should  provide  improved  ^esjlts  'f  r't,  is  contaminated 
by  random  noise. 

The  benefit  of  the  moment  method  approach  over  the  continuation  method  is 
that  no  minimization  is  reouired.  It  is  only  necessary  to  repeatedly  solve  a 
system  of  linear  equations,  which  should  De  mucn  less  time  consuming. 

It  is  important  to  prove  tne  reliability  of  eacn  of  these  methods  under  two 
important  conditions:  in  the  presence  o*  random  oo’se,  ard  with  the  number  of 
modes  present  underestimated.  If  the  rout'nes  work  well  when  the  number  of 
modes  is  underestimated,  it  is  possible  to  begin  by  assigning  very  few  modes, 
and  then  slowly  increasing  tne  number  wn;'e  using  the  previous  results  as 
initial  guesses. 

To  test  the  performance  of  the  methods,  assume  an  r',t,  composed  of  three 


modes  plus  a  OC  level 

r(t)  =  0.1  *  1  e'C‘2601tcos(2.906t  *  1  - 

*  3.3  e*°-4684tcos:9.060t  *  3 


).5  e‘°-3808tcps'6.007t  *  2. 


The  natural  frequencies  used  in  constructing  r(t  are  just  the  first  three 
of  the  thin  cylinder  target,  and  the  amplitudes  have  been  chosen  to  accentuate 
the  lower  frequency  modes.  This  response  is  then  sampled  at  500  equally 
spaced  points  between  t  -  0  and  t  *  10.  Tables  1  and  2  show  the  results  of 
using  the  described  methods  in  the  presence  of  various  amounts  of  random 
noise  and  with  tne  number  of  modes  present  underestimated .  Shown  *or 
comparison  are  tne  results  trrj m  using  tne -continuation  metnod  and  stra-ght- 
forward  Prony's  metnod.  ne^e  nc'se  nas  been  created  0/  perturbing  eacn  sampled 
point  in  the  response  by  a  random  amount  not  e«ceed’ng  a  certa'n  percentage 
of  the  maximum  'absolute  va'ue  of  tne  '-esoonse  S'nce  tne  results  *3 r  the 
direct  late- time  mimm^at’on  method  and  t“e  5  -pu'  s<»  method  are  near1  v  'lert'ca' 


those  eor  tne  E-p-u’se  teci'"  ;ue  -av< 


;*ed  (r  'jr~  t "e  tat’es  f'r  i  ■  t  < 
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Table  1.  Performance  of  various  na 


It  is  obvious  that  each  new  method  performs  quite  will  both  in  the  presence 
of  random  noise  and  when  the  number  of  modes  present  is  underestimated.  As 
expected,  Prony’s  method  fails  under  both  circumstances. 

We  will  apply  these  three  new  methods  to  the  measured  pulse  responses  of 
various  complex  targets  in  the  future.  If  any  of  these  new  methods  is  proved 
to  be  more  accurate  than  the  Continuation  Method,  it  will  be  accepted  as  our 
new  tool  for  synthesizing  the  E-pulses  of  complex  targets. 


Future  3lans 


6. 

'he  following  too’cs  --ece'  <e  ~a;c*-  rt?rv  ^  ■  "e 

( 1 )  Synthesis  of  ii-ou’ses  and  single- node  e«V'3.' •  .  .i-  .  .  ~ 

targets,  including  various  types  of  v a*  t$  and  5:a_e  .er>  '-s.  *•  ■■■* 
measured  pulse  responses  o*  tneir  sca'e  node' 

We  will  apply  our  developed  tec^r-que  *qr  s /nV'es  • :  ■ n ;  *.*<?  •  '  :,e<  s'  : 

single-mode  extraction  signals  to  v-a^ous  : one ' e ■  targets  to  esfa?'-s>' 
versatility  and  accuracy,  "he  complexity  of  experiments '  sca'e  “ode's  :* 
aircrafts  will  be  increased  by  adding  fine  structures  to  the  node's  «ner 
this  happens,  the  number  of  natural  frequencies  of  the  mode'  w'''  be 
increased,  and  consequently,  the  synthesis  of  the  discriminant  s'gna’s  wi '  ' 
become  more  difficult.  We  will  test  whether  our  technique  can  nandle  t^s 
problem. 

(2)  Refinement  of  the  technique  based  on  frequency-domain  analysis. 

In  the  synthesis  of  a  single-mode  extraction  signal  for  a  target  the 
signal  waveform  and  its  frequency  spectrum  are  found  to  be  dependent  on 
the  choice  of  basis  functions.  From  the  frequency  domain  viewpoint,  it 
is  possible  to  select  appropriate  basis  functions  which  may  lead  to  an 
effective  extraction  (with  a  maximum  power)  of  the  desired  natural  mode, 
and  at  the  same  time,  eliminating  all  other  natural  modes.  We  plan  to 
improve  the  quality  of  the  single-mode  extraction  signal  using  this  concept. 
This  is  the  continuation  of  the  study  described  in  Section  4. 

(3)  Improvement  and  modification  of  experimental  systems. 

The  short  term  plan  for  the  improvement  of  our  experimental  system  is 
to  add  a  new  sampling  scope,  a  Tektronix  7854  digital  waveform  processing 
oscilloscope  (DWPO),  and  the  modification  of  associated  instrumentation.  The 
long  term  plan  for  our  experimental  system  is  to  set  up  a  free-space 
scattering  range  inside  an  anechoic  chamber.  A  ground  plane  will  not  be  used 
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3j’  new  t'-ansm  ttmg  and  receiving  antennas  need  to  be  purchased.  Some  other 
.proponents  o*  tne  system  -nay  need  to  be  altered.  To  accomplish  this  long 
ter~  o’an,  e<tra  funding  -nay  be  needed  in  the  future. 

■i  Measurement  of  target  natural  frequencies  using  current  and  charge  probes 
'ocated  on  the  target  surface. 

we  nave  measured  the  induced  current  on  a  wire  target  excited  by  a 
Russian  pulse  with  a  current  probe.  From  the  late-time  part  of  the  induced 
transient  current  it  was  possible  to  extract  several  of  the  target's  natural 
frequencies.  The  attempt  to  measure  the  induced  transient  charge  on  the 
surface  of  a  complex  target  with  a  charge  probe  was  not  very  successful.  We 
p’an  to  pursue  this  topic  in  the  future. 

5  hew  methods  for  extracting  the  natural  frequencies  of  a  complex  target  from 
its  measured  pulse  response. 

We  will  continue  to  develop  new  methods  for  extracting  the  natural 
frequencies  of  a  complex  target  from  its  measured  pulse  response.  The  new 
methods  will  include  the  "second  generation"  continuation  method,  the  E-pulse 
method,  the  discrete  late-time  minimization  method,  and  the  moment  method  as 
described  in  Section  5.  We  will  compare  the  perforation  of  each  method  and 
finally  select  the  optimal  one  for  our  application. 

(6)  New  topics  relating  to  radar  target  detection. 

One  of  the  important  new  topics  which  may  have  great  value  in  radar 
target  detection  is  the  utilization  of  the  early-time  response  of  the  radar 
return.  The  early-time  response  contains  much  more  energy  than  the  late¬ 
time  response  does.  However,  the  former  does  not  contain  aspect-independent 
parameters  like  the  latter  does,  ’’herefore,  it  may  be  extremely  difficult  to 
develop  an  aspect-independent  scheme  using  the  early-time  response  of  the 
target.  We  may  consider  this  problem  in  the  future  as  a  ’onq  term  project. 
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Abstract— .Kb  iirki  tadcpndcal  rvdu  OffM  dbcrialudoa  icWt 
based  oa  the  uiarai  Irt^MKKi  at  the  target  It  cowMered.  Aa 
nifartlM-pdH  vtvtfora  apoa  nrtcaUoa  at  a  partlnlv  coadactlaf 
larfat  retails  la  (Be  eUaiaatioa  at  sped  Tied  natural  modal  coaleal  at  tbt 
Ka tiered  field.  EAciUiioa  ot  a  dissiaulae  larfet  produces  t  noticeably 
differcat  lale-llate  response  Cooslrucltoa  of  appropriate  etliactloa-paisr 
*a*cforau  it  discussed .  at  *HI  aa  (Be  effects  of  ra tsdom  noise  oa  (Bcir 
appiicatioa  to  thia  cylinder  targets.  Also  presented  la  cipertaoenuj 
renfkatioa  of  iBts  dtscriaiiaaiioa  coacept  utiag  umpiified  aircraft 
modeia. 

I  Introduction 

RXDaR  TARGET  identification  methods  using  the 
ime-domain  response  of  a  target  to  a  transient  incident 
waveform  have  generated  considerable  interest  recently  (l]-[4j. 
One  of  the  most  intriguing  schemes  involves  the  so-called 
“kill-pulse"  technique  as  first  described  by  Kennaugh  [5]  A 
Kill  pulse  (Af-pulse)  is  an  excitation  waveform  synthesized  in 
such  a  way  as  to  minimize  a  transient  scattered  field  response. 
Target  discrimination  results  from  the  unique  correspondence 
of  a  /f-pulse  to  a  particular  target;  excitation  of  a  dissimilar 
target  yields  a  "larger"  response. 

This  paper  describes  a  related  "extinction-pulse"  (£-pulse> 
concept,  based  on  the  natural  resonance  of  a  conducting  radar 
target  via  the  singularity  expansion  method  (SEM)  [6],  The 
time  domain  electric  field  scattered  by  the  target  is  divided  into 
an  early-time,  forced  response  period  when  the  excitation 
waveform  is  traversing  the  target,  and  a  late-time,  free 
oscillation  period  that  exists  after  the  excitation  waveform  has 
passed  [7].  (8).  The  early-ume  response  is  not  utilized  due  to 
its  complicated  nature.  The  late-time  response  can  be  decom¬ 
posed  into  a  sum  of  damped  sinusoids  oscillating  at  frequen¬ 
cies  determined  entirely  by  the  geometry  of  the  target  An  £- 
pulse  is  then  viewed  as  a  transient,  finite  duration  waveform 
which  annihilates  the  contribution  of  a  select  number  of  these 
natural  resonances  to  the  late-time  response  A  related  target 
identification  scheme  based  on  natural  target  resonances  has 
been  examined  by  Chen  [9j. 

Target  discrimination  using  this  SEM  viewpoint  is  easily 
visualized.  Each  target  can  be  described  by  a  set  of  natural 
frequencies  An  £-pulse  designed  to  annul  certain  natural 
resonances  of  one  target  will  excite  those  of  a  target  with 
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different  natural  frequencies,  resulting  in  a  different  scattered 
field  Also  made  apparent  is  the  aspect- independent  nature  of 
the  £-puise.  Since  the  values  of  the  target  resonance  frequen¬ 
cies  are  independent  of  the  excitation  waveform,  the  £-pulse 
will  eliminate  the  desired  natural  modal  content  of  the  iate- 
time  scattered  field  regardless  of  the  orientation  of  the  target 
with  respect  to  the  transmitting  and  receiving  antennas 

It  is  important  to  note  that  the  £-pulse  waveform  need  not  be 
transmitted  to  employ  this  concept.  It  is  assumed  that  an 
excitation  waveform  with  finite  usable  bandwidth  will  be  used 
to  excite  the  target,  resulting  in  a  measured  scattered  field  with 
the  desired  (finite)  modal  content  The  £pulse  can  then  be 
convolved  numerically  with  the  measured  target  response, 
yielding  results  analogous  to  £-puise  transmission  If  the 
maximum  modal  content  of  the  target  scattered  field  can  be 
estimated  from  the  frequency  content  of  the  excitation  pulse, 
then  the  £-pulse  waveform  can  be  constructed  to  yield  a  null 
late-time  convolved  response 

After  an  initial  presentation  of  the  SEM  representation  of 
the  backscattered  field  excited  by  a  transient  incident  wave  and 
calculation  of  the  corresponding  impulse  response,  two  ty  pes 
of  £-puises,  forced  and  natural,  will  be  discussed  The  results 
are  then  specialized  to  a  thin  cylinder  target,  which  has  an 
impulse  response  amenable  to  analytic  calculation.  Target 
discrimination  using  the  natural  thin  cyclinder  £-pulse.  as  well 
as  the  effects  of  random  noise  are  also  investigated  Lastly, 
experimental  verification  of  the  £-pulse  concept  is  presented. 

H  Backscattered  Field  Excited  By  Transient  Incident 
Wave 

A  perfectly  conducting  radar  target  is  illuminated  by  a  plane 
electromagnetic  wave  as  shown  in  Fig  1  The  electric  field 
associated  with  this  transient  wave  can  be  written  in  the 
Laplace  transform  domain  as 

£  (7.  s)  =  sV(s)e  "  c’*  (1) 

where  7 is  a  position  vector  in  a  coordinate  system  local  to  the 
target,  j"  is  a  constant  vector  specifying  the  polarization  of  the 
wave,  k  is  a  unit  vector  in  the  direction  of  propagation,  and 
eU)  represents  the  time  dependence  of  the  incident  field  The 
current  k  induced  on  the  surface  of  the  target  is  given  by  the 
solution  to  the  transform  domain  E- field  integral  equation 


i  |V 

- 1 


I  (7‘ ,  s)(f 


S‘  .  _ 

-t  kl 
c- 


£  <i/t  LI 

r.  s)  — 

J  4t  R 


=  -  ttn! 


00l8-9:bX  85  (MXM»2d$0l  (X)  1<>85  IEEE 


930 


IEEE  TRANSACTIONS  ON  ANTENNAS  and  PROPAGATION  VOL  A?  33  NO  9  SEPTEMBER  19*3 


Fig  l  [Hurruruuoo  at  4  conducting  target  t»>  *n  incident  ciectromagneuc 

wave 

where  t  is  a  unit  vector  tangent  to  the  surface  of  the  target  at 
the  point  ~on  5,  and  R  =  7  T  The  surface  current  can 
also  be  expanded  using  the  SEM  representation 

S 

R{7,  s)~  V  a,Rj7)is-s3) '  ‘  -  W i~.  s)  (3) 

a  •  i 

where  I T(r,  51  represents  any  entire  function  contribution  to 
the  current,  and  it  is  assumed  that  only  a  finite  number  of 
natural  modes  are  substantially  excited  by  the  incident  field  It 
is  further  assumed  that  the  onls  singularities  of  £(“,  5)  in  the 
finite  complex  s-plane  are  simple  poles  at  natural  frequencies 
5j  -  *  Juj,  Then.  represents  the  current  distribution 

of  the  ath  natural  mode,  and  aa  is  the  "class -1 "  coupling 
coefficient  given  by  Baum  (10] 

The  tar-zone  backscattered  electric  field  can  be  computed 
by  integrating  the  current  distribution  over  the  surface  of  the 
target  The  f  component  of  backscattered  field  can  then  be 
written  as 

e  ' 

£’. ( r,  s,  k)  =  — — —  eis)His.  k,  p  (4) 

r 

where  r  =  7  and  His.  k,  p  is  the  aspect  dependent  transfer 
function  of  the  target  Using  (3i  to  represent  the  surface 
current,  the  transfer  function  can  be  inverse  transformed  to 
determine  the  backscatter  impulse  response  of  the  target 

Because  ot  the  entire  function  contribution  to  the  cur-ent. 
the  impulse  response  will  exhibit  two  distinct  regions  The 
early-time,  forced  component  represents  the  backscattered 
field  excited  by  currents  during  the  time  when  the  impulse  is 
traversing  the  target,  it  has  a  duration  equal  to  twice  the  one¬ 
way  maximal  transit  time  of  the  target  T  The  late -time  free 
oscillation  component  is  composed  purely  of  a  sum  of  constant 
amplitude  natural  nodes  and  exists  for  all  time  f  >  IT  as 

N 

hu,  k,  f)=  V  a„(k.  fTe'1  cos  (w„r  *  oik.  p»,  i>2T 

(5) 

where  tt  has  been  assumed  that  the  entire  func'ion  makes  no 
contnbution  to  the  late  time  component  |b|  Thus  the  ' 
component  of  the  tar-zone  backscattered  field  is  given  simply 
by  the  convolution  of  the  time  domain  incident  field  and  the 
impulse  response,  and  is  also  composed  of  forced  and  freely 
oscillating  portions 

III  Isr  :i>fn"  f-PrisF  Wwffmrm  Svs'hims 

To  synthesize  an  £  pulse  for  a  panicular  target,  the 
convolutional  representation  of  the  backscattered  field  is 


written  in  integral  form  using  the  impulse  response  of  (5)  as 
EJ:a,  £)=  (  eW)hit  - 1' ,  k.  f)  dt' 


eW)  y  a,(k.  '  cos  [u,(/-r  ) 


+  0*ik.  f)l  dt  (6) 

This  response  is  valid  for  the  late-time  portion  of  the  scattered 
field,  i  >  =  T,  ♦  IT.  where  T,  is  the  duration  of  eu). 

The  excitation  waveform  becomes  an  £-pulse  when  the 
scattered  field  is  forced  to  vanish  identically  in  the  late-time 
Rewriting  (6i  and  employing  this  condition  yields  a  defining 
equation  for  the  £-pulsc 

V  ajk,  Oe'"  [,-f ,(  Tt)  cos  < *i,t±<t>nik,  P) 

*  B,(T,)  sin  i^,t~o,ik.  n»|  =  0.  t  >  Tl  =  T,  ~  IT 


where  tne  coefficients  .4  ,(7,)  and  B,{Tf)  are  given  bv 


f  IniT,)} 

!0 


cOS  *lmt 


The  linear  independence  of  the  damped  sinusoids  in  ("i 
requires  =  S,(  T.)  =  0  for  all  1  £  n  s  S 

It  is  important  to  note  that  A„iT,)  and  fl„(  T,)  are  indepen¬ 
dent  of  the  aspect  parameters  k  and  f  verifying  the  aspect 
independence  of  the  £-pulse  This  is  a  direct  consequence  of 
the  separability  of  the  terms  of  the  impulse  response 

A  physical  interpretation  of  the  £-pulse  can  be  facilitated  bv 
decomposing  the  excitation  waveform  as  shown  in  Fig  ?  as 

eit)  =  e'(p  ♦  e‘it)  (9) 

where  Orf  is  an  excitatory  component  nonvanishing  during  0 
s  t  <  Ti,  the  response  to  which  is  subsequently  extinguished 
by  er(t)  which  follows  dunng  T-  <  t  s  T,  Substituting  9  - 
into  1 8)  and  using  A, ( T,)  =  BjT ,)  =  0  yields 


Ci  f)e 


(cos  ]  dt 

(  sin  y 

.  ...  ( cos  u i„:  / 

e-  u  )e  i 

(  sin  *■„/  ) 


The  excitation  component  ol  the  £  pulse  necessarv  to  eradi 
cate  the  response  due  to  a  preselected  excitatorv  component 
can  be  constructed  as  an  expansion  over  an  appropnatclv 
chosen  set  of  linearly  independent  basis  functions  as 


e'(0-  y  c„g„it). 


Equation  1 1  Ol  then  becomes 

v  \f  _<r,xr„  =  -  r\ 


<1  s  v 


1  vj  rjyj  r;  w:  rj  w<  r'rj  ■ !?» »■  g  t  w  n»»  » ■«_» 


'.■".'•JKJWl 


ROTHWEU.  «  «/.  RADAR  TARGET  DISCRIMINATION 


fit) 


Fig  2  Decomposition  of  E- pulse  into  forcing  tnd  extinction  components. 

where 

*  T,  C  rnt  ^ 

dt' 


a jt  \  COS  Uft  ^ 

sin  u )it‘  J 

;J=  \  ‘  e*(t ■)«-'*  [C0S 

*  o  Sin  a)/f  J 


This  can  be  written  using  matrix  notation  as 


‘c, ' 

>1  ■ 

Civ 

Els 
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too 


(13a) 

(13b) 


(14) 


J — 20 


V 


9 

10 


0.0*28  ♦  JO. 9251 
0.1212  *  jl .912 
0.1491  *  J2.884 
0.17U  *  J3.874 
0.1909  ♦  j4  8S4 
0.208C  *  J5.845 
0  2240  *  J6329 
0.2383  ♦  j7. 821 
0  2522  ♦  J 8 . 80  7 
0.2648  *  J9.800 


Fi«  3 


Onenutioo  for  tiun-cy Under  dotation  and  first  ten  natural  frequen¬ 
cies 


Solving  this  equation  for  Ct,  ,  CiV  determines  the 
extinction  component  via  (11)  and  thus  the  E-pulse. 

It  is  convenient  at  this  point  to  identify  two  fundamental 
types  of  E-pulses  When  T,  >  0  the  forcing  vector  on  the  nght 
side  of  ( 14)  is  nonzero  and  a  solution  for  e'(f)  exists  for  almost 
all  choices  of  T,.  This  type  of  E-pulse  has  a  nonzero  excitatory 
component  and  is  termed  a  "forced”  E-pulse.  In  contrast, 
when  Tf  =  0  the  forcing  vector  vanishes  and  solutions  for 
e'U)  exist  only  when  the  determinant  of  the  coefficient  matrix 
vanishes,  i  t  .  when  det  [M(T,)\  =  0.  These  solutions 
correspond  to  discrete  eigenvalues  for  the  E-pulse  duration  T„ 
which  are  determined  by  rooting  the  determ. nan taJ  characteris¬ 
tic  equation  Since  there  is  no  excitatory  component,  this  type 
of  E-pulse  is  viewed  as  extinguishing  its  own  excited  field  and 
is  called  a  "natural  E-pulse 


9.0  t  I  4  0  9.0  9  0 

'♦orrr^  1 1  zed  t'«»  •  w 


10. 0  If.o 


Fig  4  Thin  cylinder  impulse  responses  for  9  *  30*  and  9  ■  60’  generated 
using  the  first  ten  natural  modes 


A  theoretical  analysis  of  a  thin  wire  target  has  been 
undertaker^  by  various  authors  [11|.  [121  Target  natural 
frequencies  are  determined  from  the  homogeneous  solutions  to 
the  integral  equation  1 2 )  The  geometry  of  the  target  and  us 
orientation  with  respect  to  the  excitation  field  are  given  in  Fig 
3  along  with  the  first  ten  natural  frequencies  The  frequencies 
are  normalized  by  re  L  where  L  is  the  length  of  the  w  ire  and  c 
is  the  speed  of  light,  and  correspond  to  a  wire  of  radius  given 
by  L  a  =  200 

The  thin  cylinder  impulse  response  can  be  calculated  by 
inverting  (4)  (12]  and  becomes  a  pure  sum  of  natural  modes  in 
the  late  time  Fig  4  shows  the  impulse  responses  of  thin 
cylinders  oriented  at  6  =  30*  and  9  -  60* .  generated  bs  using 
the  first  ten  modes  of  the  target  Note  the  distinct  early  and 


late-time  regions.  An  E-pulse  to  extinguish  the  first  ten  modes 
of  the  target  is  created  by  solving  ( 14)  using  the  corresponding 
frequencies 

Fig.  5  shows  a  natural  E-pulse  synthesized  to  kill  the  first 
ten  natural  modes  of  the  thin  wire  target,  using  a  pulse 
function  basis  set  The  duration  of  the  waveform.  T,  = 
2.0408E  c.  corresponds  to  the  first  root  of  the  resulting 
determ i nantaJ  equation  Superimposed  with  this  is  Kennaughs 
original  A-puise  The  simiiariity  is  sinking,  with  the  major 
difference  being  the  finite  duration  of  the  E-pulse  Also  shown 
in  Fig  5  is  a  forced  pulse  function  E-pulse  constructed  to 
extinguish  the  first  ten  modes  of  the  target  The  duration  has 
been  chosen  as  T,  -  2  3  and  the  excitation  component  has 
been  chosen  as  a  pulse  function  of  width  equal  to  that  of  the 
basis  functions 

Numerical  verification  of  the  thin  wire  E  pulse  is  given  in 
Fig  6  The  natural  E  pulse  waveform  ot  Fig  5  is  convolved 
with  the  30*  and  60*  impulse  responses  ot  F:g  4  and  the 
resulting  bat  kwanered  field  representations  ..re  >bvcrved  to  Se 
zero  in  the  late  time  Note  also  the  expected  nonzero  earls 
time  response  This  portion  is  jseful  sin.e  it  pros  oes  a 


i 
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°J  3  }  i  -  £  :  o  z  s  ? 

~-r"*4  '--v  •  _ 

-h> 

Fig  5  Pulse  function  based  £ -pulse  synthesized  to  ehmmate  First  ten  mcoes 
of  thin  cylinder  target  m  Natural  F  puise  compared  to  Kennaugh  s  K 
pulse  ib i  Forced  £  pulse  of  duration  T,  *  Z  3  L  c 

comparative  benchmark  assessing  the  quality  of  the  annulled 
component  of  the  response,  this  is  important  when  imperfect 
"extinction"  is  evident  'due  to  noise,  errors  in  the  natural 
frequencies,  etc  i 

The  question  of  £  pulse  waveform  uniqueness  as  the 
number  of  narural  frequencies  extinguished  V  is  taken  to 
infinitv  has  not  been  resolved  Fur  the  thin  cv Under  target 
there  appears  to  be  such  a  unique  waveform  Fvidence  is 
provided  hy  Fig  *  which  shows  natural  £  pulses  if  minimum 
duration  designed  to  extinguish  various  numbers  of  narural 
modes,  using  both  Fourier  cosine  and  pulse  function  basis 
sets  It  :s  apparent  that  these  wavetoms  are  nearlv  identical 
and  thev  appear  to  be  ^onserging  to  a  particular  shape 

The  most  ritual  parameter  neecssarv  for  f  puise  unique 


*».»  i  *  I  ?  . .»  m  n 

norwi  iicc  1 1»  t/iue ; 

(a) 


i 

•  i 

'  ’  H  •  *  C 


*9  :  it  io  -a  to  i  i 


' — a  :**. 

ib  i 

fig  *  Narurai  £  puoe  constructed  using  a  Founcr  cosm«  basis  functions 
to  ruminate  !  ^  "  M  and  iO  mvxlcs  <ti  Pulse  functn>n  basis  functions 
t  e.iminare  '  <  and  10  Twides 
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net*  is  the  finite  £-pulse  duration  T,.  This  must  converge  to  a 
distinct  value  as  S  —  ® .  For  the  pulse  function  basis  set  it  can 
be  shown  that  (he  natural  £-pulse  durations  are  given  by 


T,  *  2tN  —  ,  1,2,  •••,  1  S/sN  (15) 

where  is  the  imaginary  pan  of  the  /th  natural  frequency.  If 
this  is  wntten  for  the  thin  cylinder  as 

^-r[/-6</)]£  (16) 


where  4</ )//  decreases  asymptotically  to  zero  with  increasing 
then  the  minimum  £-pulse  duration  converges  to 


'  T,) ■ 


lun 

V  —  • 


IrN  L 
r[N-b(N) )  C 


-2 


L 

c 


(17) 


V  T  arget  Discrimination  With  £- Pulse  Waveforms 

Discrimination  between  different  thin  cylinder  targets  is 
Jemonstrated  by  convolving  the  natural  £-pulse  of  Fig.  5, 
which  has  been  constructed  to  extinguish  the  first  ten  modes  of 
a  target  of  length  L,  with  the  impulse  response  of  the  expected 
target  and  a  target  5  percent  longer.  The  result  is  shown  in  Fig. 
%  TV  late  tune  response  of  the  expected  target  has  been 
tuccess fully  annulled,  while  the  response  of  the  differing 
target  is  nonzero  over  the  same  period.  The  difference  in 
target  natural  frequencies  provides  the  basis  for  discrimination 
Vscd  on  the  comparison  of  adequately  dissimilar  late-time 
responses  of  differing  targets. 

Sensitivity  of  £-pulse  performance  (o  the  presence  of 
ancorrelated  random  noise  is  investigated  by  perturbing  each 
point  of  the  thin  cylinder  u  pulse  response  of  Fig.  4  by  a 
random  amount  not  exceeding  10  percent  of  the  maximum 
waveform  amplitude.  TV  result  is  shown  in  Fig.  9.  An 
•(tempi  is  then  made  to  extinguish  this  noisy  response  by 
oovotving  it  with  tV  natural  £-puise  of  Fig.  5.  As  expected, 
he  convolution  shown  in  Fig.  10,  does  not  exhibit  a  null  late- 
ume  response,  but  results  in  a  distribution  of  noise  about  the 
zero  line  Also  plotted  in  this  figure  is  the  convolution  of  the 
£  pulse  with  a  noisy  waveform  representing  a  target  5  percent 
•  mger  It  is  quite  easy  to  separate  the  effects  of  noise  and 
•arget  length  sensitivity,  suggesting  that  random jioise  will  not 
rterfere  with  target  discrimination. 


vt  F.<p*rimental  Verification  Of  The  £-Pulse  Concept 

Time  domain  measurements  of  complex  conducting  target 
responses  provide  the  means  for  a  practical  test  of  the  £-pulse 
onceps  TV  present  experiment  involves  measuring  the  near 
Kattered  field  response  of  a  simplified  aircraft  model  to 
’ransierx  pulse  excitation  A  Tektronix  109  pulse  generator  is 
i»ed  to  provide  i  quasirectangular  400  V  incident  pulse  of 
nanosecond  duration  Transmission  of  the  pulse  over  a  5  x  6 
ti  'inducting  ground  screen  is  accomplished  by  using  an 
magedi  btcomcaJ  antenna  of  axial  height  2.5  m,  half-angle  of 


Fig.  8.  Convolution  of  10  mode  natural  thin  cylinder  f-pulse  with  60*  dun 
cylinder  impulse  response  and  60*  response  of  a  cylinder  3  percent  longer 


Fig.  9,  Thin  cylinder  60*  impulse  response  generated  from  first  ten  natural 
modes,  with  10  percent  random  noise  added 


8*  and  characteristic  impedance  of  160  Q.  while  reception  is 
implemented  using  a  short  monopole  £-ficld  probe  of  length 
1.6  cm.  Although  the  receiving  probe  is  not  positioned  in  the 
far  field  region  of  the  scatterer,  the  resulting  measurements 
have  the  desired  modal  content  in  the  late -time  period, 
providing  the  small  probe  purely  differentiates  the  waveform 
Lastly,  discrete  sampling  of  the  time  domain  waveform  is 
accomplished  by  using  a  Tektronix  sampling  oscilloscope  (S2 
sampling  heads,  75  ps  risetime)  coupled  io  a  Radio  Shack 
model  III  microcomputer 
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Fig  10  ConvoUiuon  >1  nature  ter.  mode  t  pul«r  ^ * s  en 

impulse  peipo«*e 


Mr— gru  -npuMr  rf  «  McDummu  Dougi»  F  1 1  urc'Bft  model 

n»d  g>»e  iomiiSMM  — ur*j  'requeue— 


UvV 


i  fliJ 


y  ^  > 


F  puikt  -  <nsirvM. ted  u>  riimineie  ‘fte  %r  *rn  JofninAiti  mode* 
n  tie  >r  measured  mpene 


Fig  11  Measured  'espnese  W  <  B**<ng  ‘•J"  sir  »*f*  -'VBir-  end  ^er 
JnmiaaiM  nafuraj  'requrtv  <» 


In  this  experiment.  an  artempt  is  made  to  hscnminate 
between  two  aircraft  models  by  employing  the  E  pulse 
technique  Figs  1 1  and  12  show  the  measured  pulse  responses 
of  simplified  Boeing  XT’  and  McDonnell  Douglas  F  iH  aircraft 
models,  respectively  Each  model  ;s  constructed  it  aluminum 
and  has  a  geometry  as  indicated  in  the  figures  4, so  shown  n 


the  figure-  are  he  dominant  natural  frequencies  extracted  from 
'he  ale  'imr  portion  if  'he  response  using  a  nonlinear  least 
squares  urse  frting  technique  '  1  '|  E  pulse  wasetorms  can 
hen  br  onstrai  ted  >o  annul  'hese  frequencies 
P-jise  fine t ion  based  natural  E  pulses  cisnslructed  to  annul 
*ach  >1  'he  two  'arget  responses  are  shown  m  Figs  I  3  and  14 
D'Sc  nminaiion  between  'he  targets  is  accomplished  by  con 
•  •■is  ng  'hese  *aseforms  with  the  measured  responses  of  Figs 
and  2  F  g  *  shows  'he  convolution  of  the  F  IK  £  pulse 
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F>«  14  Natural  £-pulsc  constructed  10  eliminate  the  five  dominant  inodes 
in  the  F-18  measured  response. 


time  in  nanoseconds 

Fig.  16.  Convolution  of  the  707  £- pulse  with  the  F-18  measured  response 


tine  in  nanoseconds 

I).  Convolution  of  the  F-18  f-pulae  with  the  F-18  measured  response 
showing  "emnguished"  late-tune  region. 


Fig.  17.  Convolution  of  the  707  f-puise  with  the  707  measured  respoav 
showing  “ extinguished' "  late-time  region. 


with  the  F-18  measured  response.  Compared  to  early  time,  the 
late-time  region  has  been  effectively  annulled.  In  contrast. 
Fig.  16  shows  the  convolution  of  the  707  £-pulse  with  the  F- 
18  measured  response.  The  result  is  a  relatively  larger  late- 
tune  amplitude  Similarly.  Fig.  17  displays  the  convolution  of 
the  707  E-pulse  with  the  707  measured  response.  Again,  the 
late-time  region  of  the  convolution  exhibits  small  amplitude. 
Lastly.  Fig.  18  shows  the  convolution  of  the  F-18  £-pulse  and 
the  measured  response  of  the  707  model.  As  before,  the 
"wrong”  target  is  exposed  by  its  larger  late-time  convolution 
response. 


VH.  Summaky  And  Conclusion 
Radar  target  discrimination  based  on  the  natural  frequencies 
of  a  conducting  target  has  been  investigated.  The  response  of 
targets  to  a  particular  class  of  waveforms  known  as  "£■ 
pulses"  has  been  demonstrated  to  provide  an  effective  mean:, 
for  implementing  a  discrimination  process  in  the  presence  of 
random  noise. 

Two  types  of  £-pulses  have  been  identified,  natural  and 
forced.  Discrimination  based  on  natural  £-pulses  and  the 
response  of  a  thin  cylinder  target  has  been  demonstrated 
theoretically. 
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using  the  singularity  expansion  technique,"  IEEE  Trans.  Antennas 
Pro  petal ,  vol.  AP-21.  no  I,  pp.  53-62,  Jan  1973. 

[13]  B.  Drachmae  and  E.  Roth  well.  “A  confinuanoa  method  for  ideocftcn- 
uon  of  the  natural  frequencies  of  an  object  using  a  measured  response." 
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Apr  1985 


t^m*  'n  nanosecond* 

Fig  18  Convolution  of  the  F- 18  £-pulse  with  the  707  measured  response. 
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Most  encouraging  are  the  experimental  results  which  reveal 
that  two  quite  complicated,  similar  sized  targets  can  be 
adequately  and  convincingly  discriminated  using  natural  E- 
pulse  waveforms.  Further  experimentation  using  more  accu¬ 
rate  aircraft  models  is  currently  being  undertaken. 
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A  Cootiauatioa  Method  for  Identification  of  the 
Natural  Frequencies  of  an  Object  Using  a  Measured 
Response 

BYRON  DRACHMAN  AND  ED  ROTHWELL,  STUDENT  MEMBER.  IEEE 


Abstract— The  MnlldcMlM  of  lb*  aalerai  fre^Mactat  of  u  object 
M>|  wimU  data  I*  as  id-coadittooed  probieia.  A  Method  tad 
•ttarUba  la  aaiv*  lb*  proMea  booed  aa  reg  elarUattoe  by  a  coatiaaailoa 
awl  bob  la  prowled.  TW  aigeritba  la  aeptted  ta  lb*  ateaaarad  respoeae  *1 
a  aabal  aircraft,  tad  lb*  aapertortty  of  tbia  aaatbed  la  Prwoy's  art  bob  ia 
lb*  prweaca  af  beta#  la  d*w**trat*d. 


I  INTRODUCTION 

The  singularity  expansion  method  (SEM)  (2)  advocates  repre¬ 
senting  the  late-time  electromagnetic  field  scattered  from  a  finite 
sized  conducting  body  as  a  sum  of  damped  sinusoids  (natural 
modes)  An  important  problem  is  to  extract  the  natural  frequencies 
s  »  a  +  /ui  from  a  time  domain  measurement  of  such  a  scattered 
field  We  desire  to  find  a  function  of  the  form 

v 

F[A , .  .  w.v.  K.  i)  *  A  *  2  tine0n'  cos  <u>nr  +  w„)  ( 1 ) 

>i*  i 

which  "best  fits"  in  some  way  the  measured  data.  This  problem 
u  known  to  be  ill-conditioned  f  1 1  j  and  the  straightforward  use 
of  nonlinear  least  squares  or  Prony’s  method  is  unreliable  We  pre¬ 
sent  a  method  of  solving  this  problem  via  a  regularization  process- 
changing  the  ill-conditioned  problem  into  a  well-conditioned  one 
with  a  sufficiently  similar  solution 

Our  goal  is  to  give  an  overview  of  the  mathematical  basis  for 
this  method  and  a  usable,  easily  understood  algorithm  for  its  im¬ 
plementation  More  detailed  analyses  and  more  elaborate  algo 
nthms  can  be  found  in  the  cited  references. 

II  ILL-CONDITIONED  PROBLEMS 

Consider  a  problem  and  an  algorithm  to, solve  the  problem 
When  data  are  applied  to  the  algorithm  a  solution  results  If  a 
small  change  in  the  dau  leads  to  a  relatively  large  change  in  the 
computed  solution,  the  process  of  obtaining  the  solution  is  termed 
"ill-conditioned  ”  In  matrix  theory,  the  conditioning  leading  to 
the  solution  of  the  equation  Ax  =»  b  is  described  by  the  condi¬ 
tion  number  of  the  matrix  A.  which  can  be  calculated  directly 
using  singular  value  decomposition  (SVD).  or  estimated  [5J 

It  is  important  to  distinguish  between  an  ill-conditioned  a  I 
gorithm  and  an  ill-conditioned  problem  If  the  problem  is  ill- 
conditioned  then  no  algorithm  will  work  well  enough  to  pro¬ 
duce  accurate  results  (Rice  |11.  p  111))  The  extraction  ot 
natural  frequencies  from  a  measured  target  response  is  such  a 
problem 

"Regularization"  is  a  technique  for  solving  „n  ill-conditioned 
problem  by  transforming  it  into  a  related  well-conditioned  prob¬ 
lem  with  a  solution  that  is  a  good  approximation  to  the  elusive 
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solution  of  the  original  problem  See  Phillips  (9] .  Tikhonov  1 1 5 1 . 
Rice  [  1 1  ]  and  Nashed  (7J  for  more  rigorous  definitions  and  dis¬ 
cussions  Our  choice  for  a  regularization  scheme  is  to  use  a  con¬ 
tinuation  method 

III.  THE  CONTINUATION  METHOD 

Let  P(x )  be  a  nonlinear  vector  valued  function  of  the  vector 
x  Solving  the  problem  P(x)  =  0  by  Newton’s  method  or  a  quasi- 
Newton  method  requires  a  good  initial  guess  for  convergence 
[13]  The  conventional  continuation  method  may  be  viewed 
as  a  numerical  technique  to  overcome  this  difficulty  Instead 
of  solving  P(x)  a  0,  construct  a  family  of  functions 

Gt(x)  =  tF{x)  +  (\  -t)(x-Z°),  0<t<\  ( 2 ) 

where  x°  is  an  initial  guess  for  the  solution  to  P(x)  —  0.  and  solve 
the  series  of  problems  r  =  0.  r , .  ry  =  1  In  an  iterative  pro¬ 
cedure,  r  =  r,  is  replaced  by  r#* ,  =  r,  +  Ar,  and  the  problem 
GT|„,(jc)  =  0  is  solved  for  x  =  5'*1  using  x‘  as  an  initial  guess 
(where  it  is  assumed  that  Ar,  is  small  enough  to  insure  con¬ 
vergence).  The  process  is  begun  with  r  =_0.  x  =  1°  and  proceeds 
to  r  *  t f  -  1  and  the  desired  solution  to  F(x)  -  0. 

If  the  original  problem  F{x)  ~  0  is  ill-conditioned,  we  choose 
to  view  the  continuation  method  as  a  regularization  procedure 
The  problem  x  -  x°  *  0  determined  by  r  =  0  is  well-conditioned, 
while  r  *  1  corresponds  to  the  original  ill-conditioned  problem. 
Regularization  is  then  accomplished  by  making  r  as  close  to  one 
as  possible  to  adequately  approximate  the  solution  to  P[x)  =  0 
while  keeping  enough  of  the  well-conditioned  term  (x  -  x° )  to 
keep  the  combined  problem  well-conditioned 

Theorems  concerning  sufficient  conditions  for  the  convergence 
of  the  continuation  method  are  discussed  in  [8]  and  [4]  AUgower 
and  Georg  [1]  also  give  a  survey  and  history  of  the  continuation 
method. 


Tt 


(_ 

Fig  1 .  Typical  path  ii  U,  r)  space 

is  parameterized  by  arc  length  s  and  the  chain  rule  is  applied  to 
(5),  we  have 
d  <  „ 

—  H(x{s).  r(s)) 
ds 

.  dx  dH  dr 

=  (grad;//)--  +  —  -  =  0.  (6) 

ds  dr  ds 


This  is  viewed  as  a  differential  equation 


where  P  is  a  normalized  solution  to  the  homogeneous  equation 


rv  IDENTIFICATION  OF  NATURAL  FREQUENCIES  BY  A 
CONTINUATION  METHOD 

We  wish  to  obtain  a  best  fit  to  the  sampled  late-time  target  re¬ 
sponse  {rt  =  r(t,)}  by  minimizing 

-r.l1  (3) 

I 

where  P(x,  t)  is  the  fitting  function  given  by  ( I )  and  x  -  (A , ,  o, . 

.  ui/v-  K)t  is  a  vector  containing  the  unknown  amplitudes, 
natural  frequencies,  and  phases  of  the  natural  modes,  and  the  dc 
level  Consider  f(x)  as  a  vector  valued  function  with  ith  compo¬ 
nent  Ffx.  r, )  and  R  a  column  vector  with  ith  component^,  Mini¬ 
mizing  (3)  then  corresponds  to  minimizing  II  f(x)  -  R  ||2  (L2 
norm).  This  problem  is  ill-conditioned  but  can  be  regularized 
by  applying  the  continuation  method  and  minimizing 

rll/(5 )  -  R  II2  -Ml  -  Ml?  (4) 

where  x°  is  an  initial  guess  It  is  assumed  that  x  and  5°  are  ap¬ 
propriately  normalized  dimensionless  quantities  With  r  fixed, 
differentiation  of  (4)  with  respect  to  the  variables  in  x  yields  the 
normal  equation 

//(.t.  r)  =  r(grad;f )(/-  R  >  M  1  -  r>(  v  -  x°)  =  0  (5) 

where  grad;/  is  the  transpose  of  the  iacobian  matrix  o  if 

We  assume  that  15)  determines  a  simple  path  in  (.?.  rl  space 
leading  from  (5°.  r0)  to  (.?'.  ry),  as  shown  in  Fig  1  If  the  path 


then  (7)  can  be  solved  using  a  standard  ODE  follower  We  can 
also  construct  a  "predictor-corrector”  follower  to  step  along  the 
path  until  the  final  point  {xf ,  ty)  is  reached,  either  at  rf  =  I  or 
when  the  problem  becomes  too  ill-conditioned  The  predictor 
step  (Euler)  isto  solve  (8)  for  (Ax,  Ar)  with  a  prechosen  value  of 
As  The  sign  of  del  (A(s))  is  used  to  choose  the  correct  sign  of  P 
that  is.  to  continue  moving  in  the  same  direction  along  the  curve 
If  .4(s>  becomes  singular  (indicating  self-crossings  or  bifurcation^ 
of  the  path)  see  Keller  [4]  for  more  advanced  techniques  Also 
Allgauer  and  Georg  1 1 1  give  details  for  incorporating  variable  step 
size  As. 

Correction  back  to  the  curve  uses  Newton's  method  with  f  he: 
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oatlaa ta  1 from  rrt 


ctUultti  ,A°.nl,Ka  u»Ln4  linear  leapt  aquaroa 


chooae  ae  I°i“<  A^.of  ,Ka 


t.-o 


DO  until  CN  *  condition  no.  is  too  lsrc#  OP  t *1 
j  (Predictor  Stspi  Eulsr ) 
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" — — • - 

NO 

I  v, 


nornallaed  vslus  of  ^ 


I - - 

'  I°i«<°+Nr  prsdlctad  vslus  of  1 
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k  i  »0 


DO 


until  1  convergea  OR  CN  too  lore* 


C N i  «C  ond ( grsdj 3 ) 


sol  vs  {grad5H)(l,'-ikn)  «  H(Ik)  for 
jk»l  with  T  h#ld  constant 


kt-k-H 


ii-1+1 


*fi. 


-.0 

»x 


contains  idtntifisd  natural  frtq u« net 


Fig.  2.  Flowchart  for  cootiiHiatioB  method. 


•+ 

constant  The  condition  number  of  (gradj/f)  is  computed  using 
SVD  at  each  step  to  test  for  termination  of  the  algorithm.  (We 
used  a  CDC  750  which  carries  13  digits,  assumed  three  digits  of 
accuracy  in  the  measured  data,  and  terminated  when  the  condi¬ 
tion  number  exceeded  1010.)  Allgower  and  Georg  {1]  contains  an 
alternate  method  for  correction. 

At  termination  the  identified  natural  frequencies  are  contained 
in  x{  A  flowchart  for  our  algorithm  is  shown  in  Fig.  2. 

V.  EXPERIMENTAL  AND  NUMERICAL  RESULTS 

The  first  example  is  shown  in  Fig.  3  and  demonstrates  the  ad¬ 
vantage  of  using  the  continuation  method  over  Prony’s  method 
m  the  presence  of  random  noise.  The  theoretical  impulse  response 
of  a  thin  wire  at  a  30°  aspect  angle  has  been  computed  by  SEM 
using  the  first  eight  natural  frequencies  [14] ,  (2]  This  is  shown 
in  Fig.  3(a).  Fig.  3(b)  shows  the  same  reponse  with  10  percent 
random  noise  added  (10  dB  max  signal/max  notse).  Fig.  3(c)  dis¬ 
plays  the  natural  frequencies  extracted  from  the  noisy  response 
using  the  continuation  method  and  Fig.  3(d)  by  Prony’s  method. 
Clearly  in  this  example  the  continuation  method  yields  more  ac¬ 
curate  results.  Further  refinements  of  Prony’s  method  are  avail¬ 


able  (such  as  overguessing  the  number  of  modes  and  averaging 
(6) ,  [10) ),  but  we  have  never  obtained  results  as  good  as  those 
obtained  from  the  continuation  method. 

The  second  example  uses  the  measured  pulse  response  of  a 
Boeing  707  aircraft  model  as  seen  in  Fig.  4(a).  Fig.  4(b)  shows  the 
Fourier  transform  (via  fast  Fourier  transform  (FFTV)  of  the  late¬ 
time  portion  of  the  measured  response.  The  seven  largest  peaks 
were  used  to  determine  initial  guesses  for  u>,.  1  <  i  <  7.  in  the 
continuation  method.  The  resulting  best  fit  to  the  late  time  is 
shown  in  Fig.  4(c).  The  experimental  set-up  yielding  Fig.  4(a)  is 
described  in  [16] 

VI  DISCUSSION 

The  numerical  evidence  shows  the  superiority  of  regularization 
via  a  continuation  method  over  Prony’s  method  in  the  presence 
of  random  noise.  Further  numerical  experimentation  [12]  has 
also  shown  that,  in  contrast  to  Prony’s  method,  underestimating 
the  number  of  modes  in  the  response  is  not  a  computational 
disaster  In  fact,  it  is  possible  to  underestimate  the  number  of 
modes  at  first,  and  use  those  results  as  initial  guesses  when  solving 
the  problem  with  more  modes  assumed  Other  benefits  include  di- 
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Pig  3  (a)  Theoretical  impulse  response  of  a  thin  wire  target  inclined  at  9  *  30*  constructed  using  the  first  eight  natural  frequencies 

(b)  Thin  wire  impulse  response  with  10  percent  (of  maximum  signal)  random  noise  added  <c>  Natural  frequencies  extracted  from  the 
noisy  response  using  the  continuation  method  0  exact  values,  x  extracted  values  id)  Natural  frequencies  extneted  from  the  noisy 
response  using  Prony’s  method  0  exact  values.  *1  extracted  values  using  eight  modes.  ♦  extracted  values  using  16  modes 
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C„  4  Measured  near-field  response  of  a  Boeing  707  aircraft  model  to  a  nanosecond  pulse  (t>\  Founer  transform  of 
‘r^pJtre^nT^n  dotlZTnirel  modes  (peaks).  (C)  Seven  mode  best  fit  from  conunuanon  method  to  l..«  ..me  portion  of 
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rectly  incorpontmg  a  dc  level  ind  the  tact  that  overguessing  the 
number  or  modes  present  merely  results  in  negligible  amplitude 
lor  modes  not  in  the  response 
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